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Summary. — We present the spectral analysis of 14 gamma-ray bursts (GRB) X- 
ray afterglows in order to investigate the properties of interstellar matter (ISM) along 
the line of sight of GRB. We carried out a simultaneous analysis of the NIR-optical 
and X-band for those afterglows with an optical counterpart too, in order to evaluate 
and strongly constrain the absorption effect on the spectral energy distribution due 
to dust extinction from GRB environment. We evaluated the equivalent hydrogen 
column density Nh from X-ray spectroscopy and rest frame visual extinction Ay by 
assuming different type of ISM composition and dust grain size distribution. From 
our analysis we obtained a distribution of the GRB rest frame consistent with the 
one expected if GRB were embedded in a Galactic-like molecular cloud. Moreover, 
values of the visual extinction estimated from the simultaneous analysis of NIR-to-X 
band favour an environment where small dust grain are destroyed by the interaction 
with the X-ray and UV photons from GRB. 

PACS 98.70.Rz - gamma-ray sources; gamma-ray bursts. 
PACS 98.58.Ca - Interstellar dust grains. 



1. Introduction 

It's now generally believed that long-duration GRB are associated with the collapse 
of massive stars (e.g. pQ) and so a dense and dusty environment, typical of star forming 
regions, is expected in the nearby of the bursts. Up to now there are some observations 
that are in good agreement with this molecular cloud-like scenario, like the emitting and 
absorbing features observed in some X-ray afterglows (e.g. [3]), the detection of a 
large amount of dust obtained from the high resolution spectroscopy of three optical 
afterglows 0] and the high rest frame visual extinction inferred in some other optical 
afterglows (e.g. [SI ) ■ Under this point of view, the non detection of ~ 40% of X-ray 
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Table I. - X-ray afterglows sample. a Galactic Coordinates (J2000), b Redshift, obtained through 
optical spectroscopy, c Optical Transient, d Galactic hydrogen column density along the line of 
sight of each GRB (from Dickey & Lockman maps }1U( ) 
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afterglows at optical wavelengths (dark GRB) would be due to the dust absorption. On 
the other hand, the spectral energy distribution of some GRB afterglows indicates that 
dust reddening is very slow ( [HI, EI)- Moreover, for the observed optical afterglows, on 
average, the estimated rest frame visual extinction is a factor 10-100 lower than the 
expected if an ISM with dust to gas ratio and an extinction curve similar to the Galactic 
one are assumed ([E])- Also, some high redshift burst showed a larger amount of gas than 
dust. In a recent work, Stratta et al.(||]) already pointed out that absorption properties 
derived from X-ray and optical afterglows spectral analysis may require a "non-standard" 
extinction. A more careful study is required to probe the ISM of GRB environment. 
In this work we present the results of a systematic multiwavelenght spectra analysis of 
a sample of GRB afterglows in order to put some constraints on the properties of the 
circumburst environment. 

2. X-ray Data Analysis 

We selected from the XMM-Newton and Chandra archives all the X-ray afterglows 
with an high signal-to-noise ratio (see Tab[IJ in order to perform a good spectral analysis. 
Standard data reduction in the 0.1-10.0 keV energy range was performed using SAS 6.0 
for the data of the XMM-Newton EPIC instrument and using CI AO 2.3 for the Chandra 
ACIS-S instrument. All the spectra were analyzed with Xspec 11.2.0. In order to get 
Gaussian statistics, thus so to ensure the applicability of the x 2 test to evaluate the 
goodness of our fits, the spectra were rebinned to obtain at least 20 counts per energy 
channel. For the EPIC spectra we performed simultaneous PN, MOS1 and MOS2 spectral 
fitting. 

2*1. The spectral model. - According to the standard fireball model, the spectrum of 
the X-ray afterglows emission is well described by a simple power law fx{E) ex E~@ x . In 
this analysis we adopted a spectral model that takes into account also for the photoelectric 
absorption both Galactic and extragalactic due to the metal rich material along the line 
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of sight of GRB: f x oc (E' ffx x e~ N H l x e _JV "Y where Nff 1 is the is the equivalent 

hydrogen Galactic column density fixed according to the Dickey & Lockman map ^H] 
and Nfl is the extragalactic contribution to the absorption, whose value was left free to 
vary. The redshift values are derived from optical observation when available; otherwise, 
the value z=l was adopted, being this the peak value of GRB redshift distribution 11 . 

2'2. Results. - We found an absorption larger than the Galactic value in nine cases 
(GRB000210,GRB001025, GRB020322, GRB020405, GRB020813, GRB030226, 
GRB030227, GRB031203, GRB040223) with confidence level up to 4cr. The rest frame 
equivalent hydrogen column density distribution obtained has a weighted average of 
< N§ >= (1.0 ±0.1) x 10 22 cm~ 2 , consistent with the theoretical peak expected if GRB 
are embedded in a Galactic-like molecular cloud (see FigEJ. 




N H [cm" 2 ] 

Fig. 1. - (Left) The rest frame equivalent Nh distribution obtained from the analysis in this 
work and the expected N§ al distribution for GRB occurring in Galactic-like molecular cloud 
from Reichart & Price (Right). 



3. Multiband analysis 

For the eight afterglows of the sample with an optical counterpart (see TatHJ) the 
analysis was extended to the optical and NIR band, in order to better constrain extinction 
properties through a fit of spectra as wide as possible. Photometric data have been 
taken from literature. To perform this large band analysis we applied some correction 
to data. Firstly, being F(t) oc t~ a we extrapolated magnitudes at the same time of X- 
ray observation, adopting the measured optical-NIR decay index published in literature 
and then we corrected them from Galactic extinction along GRB line of sight using the 
IRAS 100 m/i E(B-V) maps by Schlegel et al. 1998 ^21 an d by deriving the extinction 
at different wavelengths using the extinction curve parameterization taken from Cardelli 
et al. 1989 assuming Ry = AyjE(B — V) = 3.1. Finally, the magnitudes have 
been converted in fluxes using the effective wavelenghts and normalization fluxes given 
by Fukugita et al. 1995 Q3J. 
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3'1. The spectral model. - For the optical/NIR data we adopted a model composed 
by a powerlaw and an absorption component that takes into account the dust extinction: 
/ (A) = C \~ 2 ~P e(~ A ( Xr ^ Avr "> Avr , where Ay r is the visual extinction. Assuming the slow 
cooling regime as the most probable one for the electron population at the observation 
time, we assumed that the electron index p is twice the X-ray spectral index previously 
obtained within 90% confidence level. We fitted the data using both the possibility 
predicted for this regime: (3 = (p - l)/2 if v opt/NIR < v c , @ = p/2 if v opt/NIR > v c , 
where v c is the cooling frequency. We tested different dust composition and dust-to- 
gas ratio using different extinction curves, assuming that the dust grains distribution 
n(a) is described by a simple power law dn(a) oc a~ q da. The curves are: the Galactic- 
like (G) from Cardelli et al. 1989 for which q = —3.5, a m j„ = 0.005/im < a < 
a max = 0.25/im ; the Small Magellanic Cloud-like (SMC) from Pei 1992 ([TB]), that's 
the same model but with 1/8 of the solar metallicity; two extinction curves (Ql and 
Q2) obtained by Maiolino et al 2000 (|32|), from simulation from study of a sample 
of AGN, for which q — —3.5, a rn in — 0.005/im < a < a max — 10/im and q = —2.5, 
a-min = 0.005/im < a < a rnax = \\im respective. We have also tested an extinction curve 
(C) derived by Calzetti et al. 2001 ^H] for a sample of local starburst galaxies. For all 
this different ISM model the rest frame visual extinction has been estimated for all the 
afterglows sample. 

32. Results. - For GRB991216 and GRB011211 we have obtained only the Galactic 
contribution to extinction, a result consistent with the previously X-ray analysis. 
For GRB000926, GRB020405, GRB020813, GRB030226 and GRB030227 the best fit 
is obtained assuming an ISM with dust grain distribution skewed toward large grains 
(Ql curve for the first case, Q2 curve for the other four cases). We also compared 
the best fit additional Nh density at GRB's redshift with the best fit Ay obtained 
for the different extinction curves. The Nh/Av relations have been compared with 
the corresponding theoretical one: Galactic, Nh/Av = 0.18 x 10 22 crrT 2 $ ]; Small 
Magellanic Cloud, N H /A V = 1.6 x 10 22 cm" 2 Q2| ; Ql, N H /A V = 0.7 x 10 22 cm" 2 and 
Q2, Nh/Av = 0.3 x 10 22 cm -2 ^7j. For the starburst galaxies no Nh/Av relationship 
has been derived due the complexity of this kind of ISM. Assuming a Galactic-like and 
SMC ISM, this ratio are well above the expected values, confirming previous studies (e.g. 
0,0)- A good agreement is obtained with an ISM with dust grain size distribution 
skewed toward large grain (see fig. Such a dust composition reconciles the typical 
low reddening observed in the optical afterglows SEDs with the high amount of dust 
observed through optical spectroscopy £T. . This kind of environment can be obtained in 
two ways: small dust grain can coagulate into larger producing a dust distribution biased 
towards larger grains as expected in high density medium (e.g.|20j). or the physical state 
of the gas and dust are modified by the intense X-ray and UV emission from GRB, (e.g 

EH)- 

Such an uncertainty will be solved by monitoring the very early afterglows and the 
possible evolution of the extinction effects. The SWIFT satellite and robotic telescopes, 
such as REM, will be very useful to provide us these informations. 

* * * 

For a more detailed analysis see Conciatore et al 2005 in prep. 
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Modello Q2 : N h z /A v = 0.3 x 1 22 cr 
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Fig. 2. - The best fit additional column density against best fit visual extinction obtained using 
Q2 curves. 
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